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Abstract. In this article we extend the results of ^] to the multivariate setting. In [T], we determined which 
families of biorthogonal functions arise as limits from the elliptic hypergeometric biorthogonal functions from 
|14l when p — > 0. Here we show that the classification of the possible limits of the BC'„ type multivariate 
biorthogonal functions from 1111 and 1101 is identical to the univariate classification. That is, for each 
univariate limit family there exists a multivariate extension, and in particular we obtain multivariate versions 
for all elements of the q-Askey scheme. For the Askey- Wilson polynomials these are the Koornwinder 
polynomials, and the multivariate versions of the Pastro polynomials form a two-parameter family which 
include the Macdonald polynomials. 



In [T] the authors considered the hmits as p — >■ of the univariate ehiptic hypergeometric biorthogonal 
functions studied by Spiridonov and Zhedanov in |14j , |12| and |13| . It turned out there are 38 distinct famihes 
of biorthogonal basic hypergeometric rational functions which appear as limits. Moreover, the degeneration 
scheme included as a subset the g- Askey scheme of orthogonal basic hypergeometric polynomials. 

In [To] and [11] the second author considered a i3C„-symmetric multivariate extension of these biorthogo- 
nal functions and established their basic properties. These properties include generalizations of all Macdonald 
conjectures: Explicit formulas are given for the squared norms, evaluation of the biorthogonal functions at 
suitable geometric sequences, and evaluation symmetry between the spectral and geometric parameters. 

In this paper we want to determine what the possible (basic hypergeometric) limits of these multivariate 
biorthogonal functions are as p — ^ 0. These limits of course depend on how the parameters (other than p 
itself) behave as p — ^ 0, so to make the question more explicit we impose conditions on how they depend 
on p, similar as in ilj. Indeed, if we write TZ^^\zi;tr;Ur;q,t;p) for the biorthogonal functions (as in jllj). 
where we have n variables Zi, four parameters tr and two parameters Ur, we set Zi — ;> Zip'^, tr — >■ trp"'' and 
Ur — >■ Urp'^'' , while forcing q and t to be independent of p. It should be noted that all variables Zi exhibit 
identical behavior as p — >■ (i.e., there is just one The limit we arrive at will now obviously depend on 
the values of a.r and 7r-. 

As in |T] there are essentially two conditions on the limits which determine whether the limits themselves 
also form a family of biorthogonal functions. First of all we must ensure that the limits are still 2;-dependent 
(and form a linearly independent set of functions), and secondly we must be able to take the limit in 
the squared norm formula. The major part of this paper consists in showing that the condition that the 
limits are ^-dependent is equivalent to the condition that the related univariate biorthogonal functions are 
z-dependent. It is straightforward to see that the condition that the squared norm formula has a proper 
limit is identical to that condition in the univariate case. As a corollary we obtain that the degeneration 
scheme of the multivariate biorthogonal functions is identical to the degeneration scheme in the univariate 
case derived in |lj. In particular, for each of the limits obtained in [1] there exists a multivariate analogue. 
We give explicit measures for which these families are biorthogonal in [2]. 

One might expect that there also exist interesting limits if one allows t to vary with p. We refrain from 
considering those cases here, as our choice of fixing t ensures that the combinatorics of the multivariate limits 
equals that of the univariate limits as discussed in the previous paragraph. 

It comes as no surprise that the Koornwinder polynomials appear as the multivariate analogue of the 
Askey- Wilson polynomials, and we obtain a multivariate analogue of all other families of polynomials in the 
q- Askey scheme. All of those families of multivariate orthogonal polynomials can be obtained as limits of 
the Koornwinder polynomials. Another interesting special case arises as the multivariate analogues of the 
Pastro polynomials. These are biorthogonal polynomials, but specializing them in a proper way reduces 
them to the original Macdonald polynomials. 
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In [TU] and [TT] two different approaches are used to define the multivariate biorthogonal functions. In 
they are defined using difference and integral raising operators, in jlOj they are expanded in terms of 
interpolation functions, which themselves are defined by vanishing conditions. Neither of these two methods 
is well-suited for taking limits directly. Thus the approach we take in this article is to define the interpolation 
functions using the branching rule. This gives us an explicit expression of the biorthogonal functions as a 
finite sum of products of theta functions. We want to take the limit in this expression by taking the limit of 
each summand, however this will sometimes lead to unwanted cancellation (where the sum of the limits of 
the summands vanishes) . We will show that the symmetries of the biorthogonal function allow us to find at 
least one expression for the biorthogonal function in which this does not happen, thus giving us an explicit 
method of finding the limits. 

The article is organized as follows. We start with a section on notation, followed by a section giving the 
definition of the interpolation functions and the associated generalized binomial coefficients. In Section [3] 
we give the definition of the biorthogonal functions and its most important properties. Section |4] describes 
how we can obtain limits of the interpolation functions. The next section proves the main result about the 
limits of biorthogonal functions: They exist as linearly independent functions of z only if their univariate 
counterparts do. In Section [5] we consider the specific limit to multivariate Pastro polynomials and their 
special case, the Macdonald polynomials. The appendix gives explicit expressions for the limits of the 
interpolation functions. 

1. Notation 

1.1. Univariate g-symbols. We say a function f{x;z) is written multiplicatively in x if the presence of 
multiple parameters at the place of x indicates a product; and if ± symbols in those parameters also indicate 
a product over all possible combinations of + and — signs. For example 

n 

f{xi,X2,.-.,Xn;z) = J|/(a;i;z), 

f{x'^^y^'^;z) = f{xy; z)f{x/y; z)f{y/x; z)f{l/xy; z). 
Now we define the g-symbols and their elliptic analogues as in [3]. Let < \q\, \p\ < 1 and set 



oo m— 1 



{x;q) - xq''), {x;q)„r^ '[[{l-xq''), {x;p, q) ^ {1 



r s\ 

xp q ) 



r— r,s>0 
m — 1 



e{xq'^;p), T{x;p,q)= II . 

All these functions are written multiplicatively in x. Note that the terminating product (a;; g)„ is also defined 
if kl ^ 1- Likewise 9{x] q;p)„i is defined for all q, though we must still insist on \p\ < 1. 

1.2. Partitions. We use the notations of [10] for partitions, which is the notation from Macdonald's book 
[5] with some additions. If A C m" then we write m" — A for the complementary partition, given by 



(m" - A), 

Moreover, if £(X) < n we define 771" + A to be the partition 

(m" + A) 



m — A„+i_i 1 < i < n 
i > n 



m + Ai 1 < i < n 
i> n 



Similarly, if Ai < to we define to" • A to be the partition 

m 1 < i < n. 



Stated differently: to" • A = (m" + A')'. 



Ai— n i > n. 



We define the relation by setting k -<m A if k C A C m" + k for sufficiently large n (for example 
n — max(£(K), £(A))). Likewise we define by setting k A if and only if k' -<m A'. In both cases we 
omit the subscript if m — 1. Recall that a chain = A'"-' ^' A^^' ^' A^^-* ^' • • • W A*-"' corresponds to a 
semistandard Young tableau of shape A'-"-', where the entries of the diagram are determined by writing k in 
the boxes in the strip A'-'^-' /X^''~^\ 

Some convenient numbers associated with A are 



i 
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Here we use j)gA' "which means we sum over all boxes in the Young diagram, i.e. we sum over 1 < i < 1{X) 
and all 1 < J < Ai. A similar notation is used for products. 

In the entire article we will use n for the number of variables Zi, which means that our partitions usually 
satisfy £{X) < n. (From context it should always be clear when we use n as number of variables and when 
we use it for the function n(A).) 

1.3. Multivariate ^-symbols. A meromorphic function f{zi, . . . , z„) is called a i?C„-symmetric p-abelian 
function if it satisfies 

• / is invariant under permutations of the z^; 

• / is invariant under replacing any one of the Zi by l/zf, 

• / is invariant under replacing any one of the Zi by pzi. 

We define the space A'-"-\uo;p,q) as the space of all _BC„-symmetric p-abelian functions / such that 

n ^{pQzf^/uo; q;p)mf{- .., z„ TT _ ^^2f^— /(. . . , z,, . . .) 

r=i f=i r(uog""^f ) 

is holomorphic for z £ (C*)" for sufficiently large m. That is, / can only have poles at the points Uoq~^p^ and 
Ug^q^p'^ for fc e Z and 1 < I < m, and these poles must be simple. It should be noted that A^^\uo;p,q) = 
A{uo;p,q) as defined in 1 . 

Let us now define the C-symbols (also written multiplicatively in x). 

(1) C7°(x;q,i;p)= [] Oiq^-H'-^x; p) Cl{x-q,t) = [] [l-q^-H^-^x) 

(2) C^{x-q,t-p)^ n eiq^'^-H^'-^x-p) C^ix;q,t)= {I - q^'-'t^'^-'x) 

(3) C+ix;q,t;p)= H 9{q^'+^-H'-^'^-^x;p) C+{x;q,t)^ l[ {I - q^^+^'H^-^'^-^x) 

The elliptic C\ are as in [TU], while the C\ are the C\ from [5] . Finally we define the A-symbols by 

.0. Clib;q,t;p) 



Alia I b;q,t;p) = 



Cl{pqa/b;q,t;p) 
which is written multiplicatively in b and 

C°^2{pqa]q,t]p) 



Ax{a I bi,...,br;q,t;p) = A^(a | 6i, . . . , 6^; g, p) ^_ ^ 



Ca {pq,'t;q,t;p)C^{a,pqa/t;q,t;p) 

which is emphatically not written multiplicatively. Here 2A^ denotes the partition with (2A^)i = 2(A|-i/2])- 
A g-hypergeometric version of Aa is defined by 

^ C°,.(ag;g,t)(7°(r;g,t) ( 1 V^' 3»(V).5»(a) 

^ ^ ^ CO(agA«)C',-(g,t;g,t)C+(a,a<zA;g,t) V aVt"-V 
Whenever no confusion is possible we omit the ] q,t]p or the ; q, t from the arguments. 



The Ca's are multivariate analogues of the theta Pochhammer symbols, while the C\'s are multivariate 
versions of g-Pochhammer symbols. The and Ax correspond univariately to the summands of a very well 
poised series, indeed 



1 - a {q;q)i \ a?q 

1.4. Transformations of generalized g-symbols. It is convenient to write down a few elementary trans- 
formation formulas for these functions, analogues of some identities for theta Pochhammer symbols. The fol- 
lowing expressions can all be obtained from the two elementary symmetries 9{px;p) = 9{l/x;p) = —^9{x;p). 

(4) Cl(j,x;q,t;p)^Ctil/x;l/q,l/t;p)^Cl^^,.^p{x;q,t;^ q-^i^h^W^ 

(5) (px; g, t; p) = C^{l/x- l/t; p) = C^.^^^,,^{x) [~ - j H-^^^^ , 

(6) Ctipx; q, t;p) = C+{l/x; 1/q, l/t;p) C+^^.Jx) (^-^) g-3"(A')i3n(A) ^ 
Likewise we can find shifting formulas for the A functions: 

(7) A", (a \pb,..., V,, ...)= A° (a | 6, .) (^^) q-^-(^')f-W 

(8) AO I 6i, . . . , br) = AO (a I 6i, . . . , 6.) ( 

(9) A;, (a I p6, . . . , z;,, . . .) = A;, (a \ b, . . . ,v,, . . .) ( -) ' ' q-^-i^')f-W 

(10) A,{^\b,,...X)=A,{a\b,,..., br) (^P^^Jl^^ ,(2-)n(A')i(-2)„(A) 

We'd like to remark that A" (a | 6i, . . . , br) is invariant if we multiply each individual bj by an integer multiple 
of p, while keeping the product Y[r fixed. Moreover, if r is even, then A" is invariant if we multiply a and 
the bj's by integer multiples of p, as long as the balancing condition JJ - bi — (apq)^/^ holds (both before and 
after the p-shift). Similarly, as long as the balancing condition pgHi — t{^Pq)^ holds A\{a \ bi, . . . , b2k+2) 
remains invariant under multiplication of the parameters by integer powers of p. 

1.5. Povi^er series in p. Most functions we are interested in are elements of the field M{x), defined in [TJ 
Section 2]. This is a field of (multivariate) meromorphic functions in the variables x — {xi,X2, ■ ■ ■), which 
can be expressed as power series / = X^teT ^'^^ some discrete set T, which is bounded from below, 
with coefficients at, which are rational functions in x. The valuation of such a series is val{f) = mintg^i 
and the leading coefficient is given by lc{f) = ayai{t) ■ Since we are interested in the behavior as p — > 0, we 
think of the valuation as describing the size of / as p — >■ 0, while the leading coefficient gives the limit of / 
(after proper rescaling). The conditions on the space imply that this limit is always uniform on compact sets 
outside the zero-set of some polynomial in x. Moreover, due to some extra conditions placed on the rational 
functions at we obtained the following iterated limit property [1] Proposition 2.3]. 

Proposition 1.1. Let f £ M{x), write p"x — {p^^xi,p"^X2, ■ ■ ■)• Then for small enough e > and any u 
we have 

lc{lc{f){p''x)) = Zc(/(p^«x)), val{f) + e val{lc{f){p^x)) = val{f{f^x)). 
With as a corollary the following important result on the valuation of a sum of two terms 



Corollary 1.2. Let f,gE AI{x) and define h = f + g. 
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• If val{f) < val{g), then val{h) — val{f) and lc{h) — lc{f). 

• If val{f) = val{g), and there exists a u such that for all small enough e > we have val{f{p'^"x)) < 
val{g{p'^'^x)) . Then val{h) = val{f) and lc{h) = lc{f) + lc(g). 

1.6. Limits of generalized g-symbols. Of course the q-symbols discussed before are elements of the field 
M{x), and since every function appearing in this article is built using these g-symbols, they are elements of 
M{x) as well. Let us now discuss the valuations and leading coefficients of the elliptic g-symbols. 
For ordinary theta functions we have 

vamxp"-p)) = -{a}{{a}~l)~-a{a~l)^ lc{9{xp'^ -p)) ^ {) /J 

where {a} = a — [aj denotes the fractional part of a. Note that val{6{xp°';p)) is a continuous piecewise 
linear function in a. The valuations and leading coefficients of the C-symbols are direct consequences of this. 
While a general formula for the leading coefficient is easily given, it becomes rather complex as it changes 
for the different C-symbols. Thus we refer to the shifting formulas ^ to note that it suffices to give the 
results for < a < 1. We have 

(11) val{C{{xp^) = \X\{\{a}{{a}-l)-\a{a-l)), (a G M), lc{Ci{xp^) 
where e = 0, +, or — . 

To take limits of A° it is often most convenient to express it in terms of C°, and take the limits of the 
C°'s. One of the important reasons we so often use the is that it is elliptic (under the balancing condition 
given above). After taking the limit, we cannot shift by p anymore, so ellipticity becomes a non-existent 
concept, thus diminishing the usefulness of this notation. 

As for Aa we'll only consider lS.\{ap°' \ t"; g, t]p). It turns out that every instance of Aa we encounter has 
t" as one of its 5-parameters. Moreover the quotient of any A^ and this one is a A|^ and we can express its 
limits in terms of C^'s as described above. Thus writing down the valuation and leading coefficient of this 
specific Aa suffices to be able to obtain the limits of the general case. We assume ^(A) < n, as otherwise 
A A (a I t") = identically 

val{/\x{ap°' I r) = -2a|A|, (0 < a < 1), 

rAi")(a;g,t) a = 0, 

Zc(AA(ap" I t")) = <^ C'O,,^, / ^ .\x\ ,,,,,,,,, 

I, {q,t-q,t) \ ^Wt^J ^ 

We would like to finish this section by making the following observation. Notice that the leading coefficients 
of these terms, only depend on whether a = OorO<Q!< 1. For general a it then follows that the leading 
coefficients Zc(C|(p"x)) and /c(A(ap" | t")) only depend on a through the component of R which contains a 
if we cut M at the integers (i.e., write R = Z U Unez("'i n + Moreover, the leading coefficients associated 
to two a's in different components, which are related to each other by an integer shifts (i.e. either both a's 
are integers, or both are non-integers), differ by a monomial factor (in x, q and t). 



2. Interpolation functions and binomial coefficients 

In this section we give a recursive definition of the interpolation functions, binomial coefficients and 
biorthogonal functions from [1^ and [11] . We use a recursive definition here as this treatment is more suited 
for taking limits. We omit most of the proofs; all proofs are given in the referenced articles. Some of the basic 
results follow from a simple recursion (together with often somewhat tedious calculations). Many deeper 
results, however, are quite difficult to prove in the presentation we give here, so we are happy to simply 
observe that our functions equal the functions defined in the cited papers, and thus must satisfy the same 
properties. 

We begin by defining 
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Definition 2.1. The binomial coefficient (^) vanishes unless u -<' X, in which case 

[a,t];9,t;p ^ ' 

11 QI'„ii,-i4.fJ-'i-i ^ 11 



t]:q,t;p 



I'/ [a.t];g,t;p 9{qt^^-H''-W, p) e{q'^'+^ ^ 'pqa;p) 

^ TT e{q^'-H^'-'pq-p) ^ 9iq^^+i-H'-^'-'pqa;p) 
, h d{q^^+^ a; p) , \{ 0{qt-'-H^+^'^~';p) 

And we have the important lemma 
Lemma 2.2. We have 

(12) /^\ ^(_pq„)|A|-|Ml^n(V)-„(^')i«(A<)-n(A)-|A| A\ 

[pa,t];g,t;p \/^/ [a. 

Proof. The resuh follows from the basic equation of 9{px;p) = —^9{x;p) and some combinatorial arguments 
to simplify the expression. □ 

This gives the coefficients in the branching rule for interpolation functions. 

Definition 2.3. The interpolation functions are defined recursively by setting RQ^^\~;q,t;p) — 1 and 
RT\-;q,t;p) ^OforX^O, and using the branching rule 

(13) Rl^''^^\...,Zi,...,v;a,b;q,t;p)= ^ cx,kRI''"H- ■ ■ . . .;a,b;q,t;p), 
where 

(14) Cx,. = CxAa,b,v;q,t;p) = / \ aV b I , „, ^'^-i^^ 



-/[.-,];,,;pA«(i"-if |t"a.,t"^,2f^;5,t;p)- 



This definition can be seen as a generalization of the expression of Macdonald polynomials using branching 
rules (see [5]). By expanding the sum further this defines i?^*-"'' as a sum over all chains = A*^°-' ^' 
^(1) ^' . . . ^' — which corresponds to a semistandard Young tableau. Thus this definition is also a 
direct analogue of the combinatorial definition of the Schur functions (which would correspond to the case 

CA,« = '(^'^/"')- 

The functions defined here are identical to the ones in [lOj, as follows from Corollary 4.5, Theorem 4.16 
and Definition 11 from [10 . In [3], Coskun and Gustafson define their well-poised Macdonald functions W\ 
using essentially the same formula (though it is quite non-trivial to identify the coefficients in the branching 
rule). The resulting equality is given by 

Let us list some of the basic properties of these interpolation functions (all can be easily proved inductively, 
except for the Zi Zj symmetry). 

• We have the following negate-the-parameters symmetry 

(15) R*^'^\-Zi;~a,-b;q,t;p) = R*^'^\z,;a,b;q,t;p) 

• The interpolation functions R'^^^^\zi; a,b) are _BC„-symmetric theta functions. Moreover they are 
contained in A'"'(6;p, q). 

• We have the "invert all parameters" symmetry 

(16) i?*("Hz,;i,i;i,i;p)=i?*(")(z,;a,6;<z,t;p)9-4"(^')i4"Wf 

a q t \ 

Moreover we have the following 



a2q2t2(n-i) 



Proposition 2.4. The interpolation functions satisfy the following shifting formulas : 

Rl^-\zf,pa,b;q,t;p) = (^-±-^y\-^-(^'h^-WRl^-\z,;a,b-,q,t-,p), 

Rt\z,-a,ph-q,t-p) = ^-^y q-^-i>^\^^^WR<-\z,-aMq,t-p), 
RfH^fpzf, ^a, ^b; q, t;p) - (^^i^) ' ' a, 6; q, t-p). 



As a corollary we obtain 
Corollary 2.5. Define the rescaled interpolation functions by 



(n), , , . p*(„), , r-'a Ct{'^)Cl{t-~^ab) 

{...,z,,...-a,b;q,t;p)^R, {. . . , z., . . . ; a,b)A,{-^ \ t ) ^^^^,,_,,^,^^,^,,,.^-^,y 



Then as a function of Zi we have g*^"^ e A^'^\b\q,p). Q*^"^ is elliptic in a and b (i.e. invariant under 
setting a ^ pa or b ^ pb). Moreover it satisfies the equation 

Q*\"^\- ■ ■ , Vpzi, Vpa, y/pb; q, t;p) = Q*^''\. . . , z^, . . . ]a,b;q, t;p). 

The specific normalization we chose for Q*^"'^ corresponds with the principal evaluation [10' (3.37)], which 
says 

Proposition 2.6. We have 

Ql^''\...,aq^H"-\...;a,b;q,t;p)^l. 

It should be noted that the interpolation functions are linearly independent for generic values of the 
parameters. Indeed the interpolation functions R*^^'^\zi;to,uo) with A C m" form a basis of the BCn 
symmetric theta functions / such that 

n 

Yl ^iPQzf^/uo; q;p)inf{- ■.,z^,...) 

i=l 

is holomorphic. 

Using the interpolation functions we define the generalized binomial coefficients (following [TOl Definition 

11]) 

Definition 2.7. The generalized binomial coefficients are given by 

^) =A^(^ \ t\pq,t;p)R*J-\...,V^q'H'-\...;t'~^^V^,-^;q,t;p). 

for any n > l{X),l{p). These binomial coefficients are independent of the choice of by (jl5[) . They are 
also independent of the explicit choice of n. 



(17) 



Let us give a basic symmetry of the binomial coefficients: 



P-J [l/o,l/f)];l/q,l/t;P [a,b];q,t;p 



in view of (|16p . 

Moreover the binomial coefficients are elliptic in a and b. 

Proposition 2.8. We have 



[ap,b\\q,t;p \^/ [a,b]-q,t-p \'^/ [a,bp]-q,t-p 



3. BlORTHOGONAL FUNCTIONS 

In this section we recall the definition and basic properties of the biorthogonal functions from [TO]. 

Definition 3.1. Let to, ti, t2, t^, uq, Ui, q, and t be parameters such that t^^"'^'^Hotit2t3UoUi = pq. We 
define 



(18) R-^\z^,to : ti,t2,t3;uo,ui;q,t;p) = X! ao 



AO(i"-Ho/uo I t^-Hoti,t^-Hot2,t^-^tot3,P'-HoUi;q,t;p)' 



The normalization for this definition is chosen so that the biorthogonal functions are highly invariant 
under shifts of the parameters. If we set n = 1 and A = (l) (i.e., A has one part) then this definition reduces 
to the univariate biorthogonal functions of |14j . 

Lemma 3.2. As functions of zi we have -R^"''(. . . , Zi, . . . ;to : ii, ^2, ^3! wq, wi; q, i;p) G A*^"^ (wo;p, q). More- 
over the biorthogonal functions are elliptic in the tr and u^, that is, they are invariant under multiplying 
these parameters with integer powers of p (as long as the balancing condition remains satisfied). Finally they 
satisfy the equations 

R^\\zip^^'^'ihP^'^ ■ tiP^^^^,t2P^^^\hP^^^^;uQp^^^,uip^^^;q,t;p) = R^^\zi;to : ti,t2,t3;uQ,ui;q,t;p) 
and 

(19) Rx (zt; T ■ -r^-r^ Ti — ^ — 5 tSW -^a [z,;to : ti,t2,t3;uQ,Ui;q,t;p). 

to ti t2 t3 Uq ui q t 

Proof. The biorthogonal functions are written as sums of functions in the space A^"^ (uq ',p,q), so they are in 
this space themselves. We can use and Propositions 12.41 and 12.81 to show that the individual summands 
in the definition of the biorthogonal functions satisfy the given p-shift equations. The final equation follows 
from a direct calculation using ([TO]) and ((T7]) . □ 

As suggested by their name, the biorthogonal functions satisfy a biorthogonality relation. There are two 
kinds of biorthogonality measures. For generic parameters we have a continuous biorthogonality measure, 
while, if we specialize to t^~^toti — the continuous measure reduces to a discrete one. The discrete 
version can be obtained from the continuous biorthogonality by residue calculus. The continuous version is 
[TTl Theorem 8.4 and 8.10], while the discrete version is ITT] Theorem 8.11] or [TOl Theorem 5.8] 

Theorem 3.3. For any partitions A and k of length at most n, and for generic values of the parameters 
such that t^^"'~'^Hotit2t3UoUi — pq we have 

^0 : ^1, ^2, ^3; Mo, ""i; q, t\p)iRt^\] to ■ ti,t2,t3',ui,UQ] q, t;p))toAi,t2,t3,uo,ui-q,t-p 

1 



(5a. K 



^^a(^ I tn,tn-Hot„tr^-Hot2,t"-Hot3,^,^,j^^;q,t;py 
where 

{q-qr{p-prv{t-p,qr 



{1,9) 



,^1 5*2 ,^3 ,^4,^5 :<?,i;j 



2«n! n"=i T{V-p, q) no<r<s<5 r(t"-^i.is;p, q) 



. . ^ n r(tzfz,^^p,g) " ntor(t.^f;P,g) dz, 

J[--- )9[--- ,Zi,--- ) II rf^±i.,±i.„ ^\ LL 



for parameters such that t^^'^~^H()tit2t3t4t5 = pq and functions / € A^'^'^{t4;p,q) and g £ A^'^\t5]p,q). Let 
ruf be such that f{zi)Y[i^{t4Z^^)/T{t4q''"^fz^^) is holomorphic, and define rUg likewise for g. Let — t^ 
for < r < 3 and t^ — t^q^^J and t^ = isq^™' . The contour is now taken such that C = C~^, contains 
all points of the form p^qHr (for i,j > 0) (and hence excludes their reciprocals), and contains the contours 
p'qHC forij > 0. 



^To be precise, C should be a chain representing the described homology class. 



// moreover toti = q "^t^ and thus t" ^^2^3^01*1 — pq™~^^ (implying no contour of the desired shape 
exists ) we define the inner product as 

{f,9)tQM,t2M,Uo,ni:q,t-p^ X! /(^O^""*?'"' ).9 (^0*""'?'"' ) 



A^^t2(n-1)^2 I t'-^r-Hoti,t^~Hot2,t'--Hot3,t"~HoUo, V'-HoU^.q, t; p) 
X 



AO / t"-^*! I *i pq pq pq .„ f.„\ 

ItO I to ' 110*2 ' "0*3 ' «OMl ' ^' '-''^^ 

and have the same hiorthogonality. 

Note that using A = k = we find that the inner products are normalized such that (1,1) = 1. 

The definition gives us an expansion of f6^^ in terms of the interpolation functions R*J'^\; to , uo) ■ This 
is not the only basis for A^"^{uo), indeed R'^^^\; v , uo) is such a basis for all values of v. It turns out we 
need expansions in these different bases of interpolation functions as well, as some limits of interpolation 
functions are independent of Zi (so we do not want to expand into those). It is convenient to first recall the 
definition of from lOJ. 

Definition 3.4. We define the coefficients ^Ia/k 
^x/K{a,b;vi,V2,V3,V4;q,t;p) 

^ V A\ [A ^lipqa' I W«fe)A^(f I ,^) -i^ Ctjl^Kit) 

«CMCA \pqa-^,pqah\ \^ J U I ^J^kI 6^pg I ahpq ) r=l ^M*- )^ b > 

Univariately fl\/K, corresponds to a very well poised series 12V11. In |10j it was shown that these coefficients 
satisfy a Weyl group of type 1)4 symmetry: 

Theorem 3.5. flxfi^{a,b;vi,V2,V3,V4,;q,t;p) is symmetric under permutations of vi, V2, fs and W4 and 
satisfies the equation 

^x/Kia,b;vi,V2,V3,V4;q,t;p) = flx/nia, — ^;wi,W2, 1/^3, l/vi;q,t;p) 

V3V4 

Moreover we have 

^\/K{a, b; Vi,V2, V3,V4; q, t\p) = nx/^{-a, -b; -vi, -V2, -V3, -W4; q, t;p) 

There also exists an evaluation formula for the JIa/k if the product of two w-parameters equals abpq. This 
is the bulk difference equation [lOl Theorem 4.1]. 

Theorem 3.6. We have 

O ( M a^.fA Cti'^^,pqav,,pqav2,^,,^-^)CUpVa^) 

U,/^(a,,\v,,V2,x, ^ ]-y^^^^^^^^^o(H™)co(p,az;i,p,az;2,^,f,T)- 

Finally we would like to observe that 

\a ( V1V2V3V4 I £r«2iaHi.'l 4 /^Ofv,.\ 
n,/^ ia,b;v„V2,V3,V4;q,t; p) ^0(p^^2 | ^g,^) {[ cfi^) 

is elliptic in a, b, Vr, i.e invariant under a — > ap, b — > bp, vi — > vip, etc. 

We can get other expansions for the biorthogonal functions, by using the definition and then expanding 
ii!*^"^(; io, Wo) in i?*("-)(; w, uq) (for any v of our choice) using [TOl Corollary 4.14]. This gives the equation 

(20) E}^\-to:tiMM;uo,ui-q,t-p)= ^ "° 



Cl{^„,t--Hot,,t-~Hot2,t-~Hot3,^^J 

^ r> I "OS pqs pqs pqs 

X > "A/jyls, 1 — ; r, r, T-,UQVs;q,t;p) 

^ f'-Ho totif^-^ tot2t''-^ totst""-^ 



wq^o "0*1 Uot2 upia 



where s = 1/ ^/pqu^ui. Note that the summation over v (without the terms in front) is permutation 
symmetric in toi tii ^2 and ^3 (using the symmetries of VL\/i,)\ in particular this shows that the biorthogonal 
functions themselves are symmetric under the exchange of and ti (or just any permutation of (io, tl^ t2, ts)), 
as long as we multiply by the right product of C^'s. The current normalization is such that 

(21) R'^^\r-Ho;to : h,t2,h;uo,ui) = 1. 

Together with the permutation symmetry in the tr this also gives us evaluations for i?^"'' {t"~Hr) for r = 1, 2, 3. 

For some values of v we can evaluate the factors in this expansion, not surprisingly these are v = tr 
(r = 0, 1, 2, 3), which lead to the original expansion (if v — Iq), and versions of that one using the symmetry 
in the tr- 

We would like to mention a few results on these biorthogonal functions, which are essentially algebraic 
equations in our space of formal power series in p. Thus it should be relatively straightforward to obtain 
limits of these equations. However we have not yet completed the full study of obtaining these limits. 

We begin with the following evaluation duality [TUl Theorem 5.4], which is the analog of Macdonald's 
evaluation duality conjecture. Once we have obtained the limits of the biorthogonal functions the relation 
will reduce to duality relations for our limiting biorthogonal functions. 

(22) i?l,"^(toi""'g"'; io : h,t2,h;uo,ui) ^ R^^\ke~'q^^ ; k : ii, ^2, 4; ^0, wi), 
where the new parameters are given by 

ia = V^otihts/pq, idr = toU, {r = 1, 2, 3), io/iir = to/wr, {r = 0, 1). 

Notice that we need the valuation of z to equal the valuation of to (i.e. C = Q^o) for the limit to work (though 
by permutation symmetry in the tr^s and the z ^ 1/z symmetry we have a few more choices). Also observe 
that the case k = is exactly the normalization equation. 
Let us define a difference operator [TT, (6.18)] 

^f"(";'..)/(«.):^ E %|5£f^ n |?Kf 

It is easily shown that if t'^~^voViV2V3 = p this difference operators maps the space of i3C„-symmetric abelian 
functions to itself. We can rescale this operator to 

n(„). , w D^''\vo,vi,v2,p/t'^-'^voViV2;t,p)f 
D\ >{vQ,Vi,V2;t,p)f = „ — — r 

so that we obtain the difference equation 

D':^\uoMM\t.P)Rx\-,q^'^h ■■ q^/\,q-'^^%,q~^^%;q^^^Uo,q-^/^Ui) ^ i?i"^(-;io : ^i, ^2, ^3; Mq, 
In the same vein we can define 

D~^"\uo;t,p) = D^g^\uo,quo,p/uo,l/t''-\oq;t,p) 

and we set 

n+Mf ^ \ff \ TT 0{pqt''~'vi/vo;p) 

Dt^ Hvo : V, : .3, .4; t,rt/(..) = H U2^r<.0ivrtn-^v,;p) 



where V5 is determined by the equation Y[r=o ~ P^l- Acting on the biorthogonal functions these 
operators give the equations 

Z?+(")(wo : to : tiMM)Rt\-^q^'^tQ ■ q^'^ti,q^'h2,q^'^t:i-q-^/^uo,q-'^'^ui) = R^;'l,,.{-,tQ : ti,t2,t3;uo,Ui) 
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Figure 1. The fundamental domain including its partition in a octahedron and two tetra- 
hedra. The points are labeled by the values of (a,/3, C)- The difference between solid lines 
and dots on the one hand, and dashed lines and open dots on the other, is that the former 
would be visible were the parallelepiped solid, while the latter would be invisible. 



and 

{uo)Ri% (•; q-'/'to : q-'^h, , q-'/%, q-'/Hs; q^^uo, q^'^m 



A+l 

4. Limits of interpolation functions and binomial coefficients 

In this section we discuss the limits of interpolation functions and binomial coefficients. This section only 
contains the methodology and a description of the results. The explicit calculations and expressions for the 
explicit limits we obtain are relegated to Appendix [X] 

First we wish to remark that the interpolation functions, and therefore also the binomial coefficients, 
are elements of the space M{zi,a,b,q,t), respectively M{a,b,q,t), as explained in Section [1] In particular 
finding their limits consists of determining the leading coefficients, and convergence is then automatically 
uniform on compacta in C* outside the zero set of a polynomial. 

Interpolation fimctions are defined using the branching rule (fT3|) . which is 

R*x"'^^\- --^Zi,.. .,v;a,b;q,t;p) = ^ CA,«i?*'"H- --^Zi,.. . ;a,b;q,t;p), 

SO in order to obtain their limit we first need to find the limit of the coefficients ca.^, which were given by 
p^ . As we already know the limits of the C°'s and Cj?'s appearing, we can write down this limit as soon 
as we know the limit of (^), ^, ^ . As (^) is given as a product of theta functions, that limit is immediate. 

It turns out that the valuation of cx,K,{vp'' ap" , bp^) is always of the form a;(|A| — |k|), where x = x{C,, a, (3) 
is an explicit piecewise linear function. In particular, an immediate induction shows that the valuation of 
R*x^"'\zip'';ap°',bp^) equals a:|A| (assuming no cancellation occurs), and the limit of i?^^"^ can be defined 
recursively using a branching rule (|13l) . where we use the limits of ca,k instead of ca.k themselves. 

So far, we have discussed the limit of R*^"'\zip^]ap'^ ^bp^) for a given (q!,/3, C) G M?. However, as there 
are infinitely many such vectors, we now want to bring some order in these different limits, and show that 
in fact there are only finitely many essentially different limits. 

Recall the p-shifts in the arguments of i?^*-"^ as given in Proposition 12. 4[ and the fact that i?^^"^ is 
p-abelian in the z^-parameters. As a consequence we can rewrite R*^^"'\zip'^; ap°' ,bp^) for any (a,/3, C), 
as an explicit power of a, &, q, t and p times an interpolation function with parameters {a, /3, C) in some 
fundamental domain of the group G of shifts along the lattice generated by (1,0,0), (0, 1,0) and (^, ^, \)- 
In particular, up to a simple multiplicative factor, all the limits of the interpolation functions, are limits of 



interpolation functions with (a, /3, Q in such a fundamental domain. Let us fix this fundamental domain as 
the parallelepiped 0<a±(,f3 + (< 1. 

Now recall that the leading coefficient of C^{p^x) only depends on x via the component containing x in 
the decomposition M = Z U Unezl"-' + same property holds (by inspection) for Now we 

immediately see that the leading coefficient of ca,k only depends on where the numbers (a — /3, a + C, ^ 
C, + /? — C) are with respect to this partition of M. In particular, if we cut by the hyperplanes a — (3 e Z, 
a ± C g Z and /3 ± C G the leading coefficient of R*^"'\zip'^ ; ap°', bp^) only depends on where the vector 
(a, /?, is located with respect to these hyperplanes. 

These hyperplanes cut space in a tessellation of square pyramids and tetrahedra, where the square pyra- 
mids pair up into octahedra, such that around each octahedron we only find tetrahedra and vice versa. 
Our fundamental domain consists of one octahedron (i.e. two square pyramids) and two of its neighboring 
tetrahedra, see Figure [T] The limit now only depends on which open polytope (i.e. vertex, edge without 
vertex, triangle/square without edges or interior of tetrahedron/pyramid) in this tessellation contains the 
vector (a,/3, C). Of course two polytopes related by a shift from the translation group G will provide the 
same limits (up to an explicit power of a, b, q and t). 

So in principle we have 1 limit associated to a vertex (the group G acts transitively on the vertices of 
this partition of M'^), 6 to an edge, 8 to a triangle, 1 to a square, 2 to a tetrahedron, and 2 to a square 
pyramid. We can use symmetries of the interpolation function we have not yet used to cut back this number 
somewhat. 

First of all the Zi — >• l/zi symmetry allows us to identify limits corresponding to vectors (a,/3, C) and 
(a,/3, — C)i limits which are related by a reflection in the plane C = 0. Secondly the {a,b;q,t) — ?> 
{l/a,l/b;l/q,l/t) symmetry (ITBl) . allows us to relate the limits at (a, /3, C) and (— a, — /3,C), i.e. limits 
related by a rotation of 180° around the i^-axis. In practice we will often prefer not to use the last symmetry, 
as it breaks the condition \q\ < 1, which is necessary in the measures (though not for defining the interpolation 
functions). 

Writing down the limits explicitly we find that the limits associated to vectors inside octahedra (i.e. either 
the square, or one of the two pyramids) are independent of Zi (and therefore not particularly suited as a 
basis of functions to expand other functions in) . However the limit associated to a vector in a tetrahedron 
becomes the Macdonald polynomials (in either Zi or 1/zi), and thus in particular an independent set of 
symmetric functions of the Zi. Now recall the iterated limit property, Proposition ll.il This implies that if 
we have a limit L associated to the polytope P, we can take further limits of L (by sending some parameters 
to or infinity) obtaining the limit associated to any polytope which contains P in its boundary. All limits 
not associated to the interior of the octahedron have a limit to the limit associated to the interior of one 
of the tetrahedra. In particular all limits not associated to the interior of the octahedron must form an 
independent set of functions of the Zi. 

While these families of limiting interpolation functions all form independent families of functions of the 
Zi, they do not all span the same space. Indeed, as mentioned the interiors of the tetrahedra correspond to 
permutation symmetric polynomials. If the vector (a, /3, C) is on a plane /3 + G Z, respectively /3 — C e Z, 
then the functions have poles at Zi G b~^q^^^, respectively Zi € bq^^-^, as can be seen by inspection. Thus 
whenever /3 + C€Zor/3 — the limits are rational functions of the Zi. Moreover we get limits which 

span spaces of i3C„-symmetric functions (i.e. invariant under Zi ^/zi) or functions which lack this kind 
of symmetry, and other distinctions. 

To obtain the limits for the binomial coefficients (^)j ^ fDefinition l2.7p we can now just plug in the 
known limits of the interpolation functions in their definition. The result can be described in much the same 
way as the limits of the interpolation functions. However, we consider the limit of (^) „ , and see 

that the limit depends only on which part of the plane IR^ contains (a, /?) if we cut the plane by the lines 
a e Z, /3 € Z and a — /3 € Z. It turns out the limits are quite a bit prettier as there are no constant factors 
in the p-shifts (Proposition of the binomial constants, which implies that the valuation of the binomials 
is always 0. 

Let us end this section by formulating the main points in a proposition 
Proposition 4.1. We have the following results about the limits of the interpolation functions. 
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(1) For each vector {a,l3,(^) G R"^ there exists a scale a;(a,/?, C) such that 

val{R*^"'\z,p^\ap°',hp^]q,t]p)) = | A|x(a, /3, C), 
lc{Rl^"'\z,p'^;ap°',bp>^;q,t;p)) -. Rl^^l i^^^^{z,; a,b; q,t) 

^*x^a 13 Q rational function of the Zi £ C* with poles in at most z^^ € hq^^-^ . 

(2) The translation group G generated by shifts in the directions (1,0,0), (0,1,0), and (i, i, i) when 
acting on the vector (a, /3, C), leaves R*x^l^ p q invariant up to an integer power of a, b, q, and t. 

(3) The function R*x^l^ p q depends on {a,/3,(^), only through the location of this vector with respect to 
the hyperplanes a — /3, a ± ^, /3 ± € Z. 

(4) // (a, /3, C) is not in the G-orbit of the interior of the octahedron (from Figure Q]) then the -R^^(2 i3 c) 
form an independent family of functions of the Zi, that is, there does not exist a non-trivial equation 
of the form 

^dx{a,b,q,t)Rl^'^lj;^^^{z,;a,b;q,t) = 0, 



A 



where the sum is finite. 
(5) For any vector (a, /3) G 




val n =0, Ic ^ 



[ap'',bpf<]-q,t-p/ \\l^/ [ap'',bpP]-q,t-p/ \l^J (a,lj):[a,b];q,t 

The translation group H generated by shifts in the directions (1,0), and (0,1) when acting on the 
vector [a, (3), leaves (^)^ ^ invariant. The limit (^)^ ^ depends on (a,/3) only through the 
location of this vector with respect to the hyperplanes a, f5, a — f5 € Z. 

5. Limits of the biorthogonal functions 

The goal of this section is to show two things. First of ah we want to show that the Hmit of the biorthogonal 
functions are z-dependent if and only if the limit of the univariate biorthogonal functions are z-dependent. 
Moreover the limit of the biorthogonal functions with A C m" form a basis of the appropriate limit space, 
when this occurs. Secondly we want to show that there always exists an expansion (of the form pop . with 
^x/iy replaced by its definition, Definition l3.4p for which we can obtain the limit of the biorthogonal functions 
by simply interchanging limit and sum. 

In order to study the limits of the biorthogonal functions, it is necessary to first consider the limits of 
i^x/ii- Indeed, interchanging limit and sum in the defining expansion ()18p of R^^^ does not always work. In 

that sum we expand in Rli^^\;to,uo) (for pL C A), and we have seen in the previous section that for 
some ways of scaling the z^, ^-nd uq with p the limit of these interpolation functions will not depend on 
z. This might mean that the limit is itself ^-independent, but it could also be that the valuation of R^^^ 
is more than the valuation of the summands in (I18p . By using equation (j20p instead we can expand R^^^^ 
in R*^^^\; V , Uq) for a f of our choice, and in particular we can choose v so that at least the interpolation 
functions have a proper limit, and their limits form an linearly independent family of z-dependent functions. 

5.1. Limits of Omega. In this subsection we show that at least one of the eight different sum expressions 
we have for f^A/K (using its symmetries) is such that we can determine its valuation as the minimum of 
the valuations of the summands (in the expansion in powers of p). In particular this implies we can always 
determine the valuation of ^x/k- 

Recall that the valuation of a sum of a finite number of terms is at least the minimum of the valuations of 
the individual summands, and equals the minimum if the minimum is attained at most once. Moreover, if 
the valuation of the sum equals the minimum of the valuations of the summands, then the leading coefficient 
of the sum equals the sum of the leading coefficients of summands with minimal valuation. On the other 
hand, if the valuation of multiple summands is the same, the leading coefficients of those summands could 
cancel each other, and we would have a priori no good way of expressing either the valuation or the leading 
coefficient of the sum in terms of leading coefficients and valuations of summands. 

13 



Let us first consider the case where we can determine the valuation and leading coefficient of ilx/n because 
there is a unique term in its sum expression with minimal valuation. 

Proposition 5.1. Define the function 

4 4 4 

/(a, ^; 71, 72, 73, 74) = g{a + P -^jr) + g{2a) - g{-a - /3) - 7,- - 2/3) + ^ 5(7^ - P) - g{a - 7^), 

r— 1 r— 1 r— 1 

where g{x) — {a;}(l — {x}) (with {x} the fractional part of x). Let a — ap" , b — bp^ and Vr ~ VrP'*'' (where 
a, b and Vr are all p-independent) . 

• If f{a, j3]^) > we have 

lc{nx/Mb-Vi,V2,V3.Vi]q,t-p)) ^lc\\\ ,^ I .l...3.4 ^ ^+/aN 1 ) 

and likewise with the valuations, i.e. 
val{nx/^{a,b;vi,V2,V3,V4;q,t;p))^val \ \ ^nr^Q ao ^ ^1.2^3^4 1 v,v l^ r+fa^ [ ,J 

• If f{a, l3;j) < we have 
lc{il\/K{a,b;vi,V2,V3,vr,q,t;p)) = M J_ 



11 r'0('P2£.'l AO ( a I /^+|- tin)2-U3f4 \ I / 



and likewise with the valuations. 
Proof. Recall that ([TT|) implies 



1 



= -|A.|(W({a} - 1) - a(a - 1)). 



(■ , 

yvl [a,b] 



Also recall from the previous section that val ( (^) f„ ) = 0. Thus for the individual summands in Definition 
of ^x/k we find 



^ I li riO(m±\r<OC£r^\ \Q(a 1 1 \ A0/' fii'2^'3«4 1 viV2V3V±\ \ III \ K i 

\r=l ^t^^ v-r l^l^y b I I ^l^l^y b-^pq I abpq > W [pqa^ ,pqab] \^/ 

In particular we see that if f{a,l3,^) ^ the valuation of the summands is |/i|-dependent. If /> it is 
minimized if is largest, thus only at = A. On the other hand if /< it is minimized \i = k. Together 
with the discussion before the proposition on how the valuations of summands correspond to valuations of 
the sum, this finishes the proof. 

We can simplify the summands aX p — \ and /i = k by using the explicit value of (^)j^ j^j.^ ^.^ from jlQl 
(4.8)]. □ 

In case the function / in the above proposition vanishes we still don't know the valuation of I^a/k or 
its leading coefficient. However we can revert to CoroUarv 11.21 of the iterated limit theorem to tell us the 
valuation and leading coefficient of J^a/k in some more situations. 

Proposition 5.2. We use the notation of the previous proposition. If f(a, /3;j) = 0, and there exists a 

) such that there exists a S > such that for any e with S > e > we have 
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/(a + exa, l3 + ex(,; 71 + nx^^ , . . . , 74 + ex^^) ^ then we have 
val{nx/Ja, b; vi,V2,V3, vr, q, t;p)) 

^'^ I li no ( P2±\r<o ( ihL\ AOi'a I 1 y/yo < ""i-"2f3-"4 1 V1V2V3V4 ^ I / I 

''^K^ 6 I pqab^'^K^ h^pg I ^bpg ) / [pqa^ ,pqab] 

for any k d fi d X. Moreover we have 



r=l ^"-"kV b > ^l-i\b I pqab^'^K^ ti^pq I a6pq V'"^/ [pqa^^ ,pqab] 



Note that the different expressions for the valuations above do not depend on /i precisely because / 
vanishes. 



Proof. This is a direct consequence of Corollary [TT2] extended to an arbitrary finite number of summands. If 
/(a + eXa, P + exb', 71 + £2^71 , . . . , 74 + ex^^) > (for all small positive e) then the summand with /i = A will 
dominate all other terms in the iterated limit, whereas if /(a + eXa, 13 + exb] 71 + tx-y^ , . . . , 74 + ex-y^) < the 
summand with = n will dominate all the others. □ 

Notice that g{x) — x{\ — is a piecewise linear of x. Now we can see that / is piecewise linear, by 
replacing all instances of g in the definition of / by x{l — x) and observing that the result vanishes. Hence, 
the points of space on which / is identically zero in a neighborhood around that point are some polytopes. 
To obtain a proper limit for on those sets we can use the W(D^ symmetry. Indeed, while the sum 

^A/ft is invariant under this symmetry the summands are not, and in particular the function / which controls 
whether we can obtain the limit is not. / being clearly permutation symmetric in the 7j we only have to 
consider the 8 cosets of S'4 in W(D/^). This gives us 8 functions /. It turns out only 5 of these functions are 
linearly independent. However, in each point of space at least one of these five (in fact we only need 4 of 
them) is non-zero at this point, or at least at some points in an arbitrary small neighborhood of this point. 
In particular we can obtain a proper limit for ri^/^ at all points in space. 



Lemma 5.3. Let f he as defined in Provosition \5.1\ Consider the four functions 

f{a,P;'y), /i2(a;/3;7) = /(a,/3-7i - 72; -71, -72,73,74), 
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/i3(a;^;7) = /(a,/5-7i - 73; -7i, 72, -73, 74), /i234(a; /3; 7) = 2^ 7r; -71, ^72, -73, -74)- 

Then for any point (a,/?; 7) at least one of these four functions is not locally constant and, hence, not locally 
zero. 

Proof. As / is piecewise linear we can differentiate almost everywhere; in particular we may assume that all 
four functions are differentiable at the point (a,/3;7). Moreover, if all four functions are locally constant, 
the derivative (in any direction) has to vanish. In particular we get 

= 4(5'(2a) - g'ia - 71) - g'{a + 71)). 
Now we can plug in g'{x) = 1 — 2{x} (which holds if a; ^ Z) and divide by 8 to get 

= -i - {2a} + {a- 71} + {a + 71}. 

As the sum of the fractional parts in this equation is an integer, this equation cannot be satisfied, so our 
assumption that all derivatives vanish is false. Therefore at least one of the four functions must not be 
locally constant. □ 
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The upshot is that for any vector (q;,/3;7) we can determine the size of J^a/k- Combining the previous 
lemmas we thus get 

Proposition 5.4. There exist piecewise linear functions /i(q;,/?;7) and f2{a, I3;"f) such that 

val{nx/Kiap'^,bp^;VrP''^;p,q)) = ^; 7)|A| + /2(a, /?; 7)|k|. 

In fact we even have exphcit, but very ugly, expressions for /i and f2', indeed, using whichever of the four 
functions in the previous lemma is not locally zero at a point, will provide us with the values of /i and f2- 
In specific cases (i.e. when we are given a, l3 and 7) we can with relative ease (by checking all four cases) 
find an expression for the valuation and leading coefficients of flx/,^. 

5.2. Limits of Biorthogonal functions. Now we consider what limits we can obtain from the biorthogonal 
functions. We consider here limits of 

(z,/ ; top"" ■.tip''\ i2P"=^ , i3P"^ ; uop'"" , uip^^ ■,q,t;p), 

where the parameters Zi, tr, Ur, q, and t are independent of p. 

If we consider the expansion (|20p in terms of we can determine the valuation of the summands, in 

particular the summands have valuation hi\X\ + h2\h'\ for some piecewise linear functions hi and /12 of 
a and 7. We can now use a similar method as for to determine the valuation of the limit. However 

because we know the limits of the interpolation functions often form a linearly independent set of functions 
of Zi we have an extra way of seeing if cancellation might occur when the valuation of all summands is equal. 
This allows us to simplify the argument. 

Proposition 5.5. Suppose f^("~"'^'ioii^2i3'"o'"i = PQ- Write t^ — trp"^'' , Ur — Urp'''' , Zi = Zip'^ , and v — vp'^ , 
where tr, Ur, Zi, t, q, and v are independent of p. Let hi and /12 be piecewise linear functions such that (with 
s = 1/y/pquQUi) 

val ilx/^ls, ; ———I, , , , ——^.uovs; q,t;p)C^{—-, — -, — -, —-)Rj^ >[zf,v,uo) 

\ t" Ho totit" ^ ^0*2* ^0*3* Uoto Uoti Uot2 Uotz 

= /ii(C;a;7;'^)|A| +h2{C,;o.\-i]v)\v\. 

• // /i2 < then 
lc{R^^\zi]to ■■ ti,t2,t:i;uo,ui)) 

(nO/pq£zll2. PI PI PI PI \ C'^ I 1 \ 

C0(^,t"-1^0^„t«-t0^2,^"-^0^3,^)C,+ (^:^)''^ 

and the valuations are also equal. 

• If ^2 = and (C;i^;7o) is not in (a shift of the) interior of the octahedron (from Proposition 
then 

lc{R^^\z,;to : tiMU.uo^ui)) = V^cf ,f ' ^ ^ uptp ^ .pu ' ^o*. ■ W 

I "0'^ P^^ P1^ P1^ + \U*in)f \\ 

X il\/^(s, T — : -, -, T-.uoVs; q, t; p)K,} '(Zi\ v, itn) 

' t''-Ho totit"^-^' tot2t'^-^ totst''^^ ' ^ 7 

and the valuation of the biorthogonal function equals the valuation of each of the summands. 

• If h2 > then 



lc{R'-^\zi;to : ii, ^2, ^3; "0, wi)) = Ic' ^ 



A V n.n > 



C'xi^,t--'toti,t^'-Hot2,t--Hot,^j^^) 

uqs pqs pqs pqs 

X S2a/o(s, — ; T-.uovs; q, t;p) 

' f^-^to totit"--^ tot2t"--^ totat"-^ 



and the valuations are also equal. 
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It should be noted that we do not cover all possible cases in the proposition above, in particular we do 
not claim anything for the case h2 — and (C; i'; 70) in (a shift of the) interior of the octahedron (from 
Proposition 14.11) . 

Proof. If /i2 < then the valuation of the summands in the expansion (PH)) of the biorthogonal function in 
terms of is minimized at v = X, so we know the leading coefficient is just the leading coefficient of the 
ly — X term. We simplified that term by using that is a sum of just a single term, and again using the 
equation for (^) . If on the other hand ft,2 > then the leading coefficient of the biorthogonal function equals 
the leading coefficient of the — summand. 

Finally if /i2 = the valuations of all summands are equal, however if (^; 70) is not in (a shift of) the 
interior of the octahedron, we see that the sum of the leading coefficients of the summands cannot vanish, 
as the limits of the interpolation functions form a linearly independent family (as functions of z). Thus the 
valuation of our biorthogonal function must equal the valuation of the individual summands, and the leading 
coefiicient equals the sum of the leading coefficients of the summands. □ 

Now note that we can always choose v such that (C;j^;7o) is not in the interior of the octahedron, for 
example by taking 1/ — In those cases this proposition exactly tells us what the limit of the biorthogonal 
function is. In particular, in that case, we see that the limits form an independent family of functions of the 
Zi if and only if /12 < 0, which is equivalent to the condition that the limit of R^^^^ depends on Zi at all (for 
any A 7^ 0). Thus we can determine which limits form such a family (and therefore have a shot at being one 
part in a pair of biorthogonal functions) by considering whether the limits of the univariate biorthogonal 
functions with those parameters depend on z. And this is precisely the situation we studied in [1] Section 
4]. 

When writing down explicit limits, it is better to first try whether the cases v = and = ar, (0 < r < 3) 
work, as this would lead to an expression of the limit as a single series (over partitions) (using the evaluation 
from Theorem 13. 6|) . It might be that the vector {C,\ar\^a) is always in the interior of the octahedron, in 
which case we have no choice but to use a different v and expand the leading coefficient as a double series, 
but fortunately this happens only very rarely. 

The final question is when the limit of these families of biorthogonal functions still form a biorthogonal 
family of functions. This can only be true if both families of biorthogonal families are still z-dependent, 
and if the valuation of the norms is correct. This condition is easily seen to be equivalent to this condition 
in the univariate case. Therefore the question when the limit can still form a biorthogonal system reduces 
to the same combinatorial problem as in the univariate case. In particular [1] Theorem 5.2] also holds for 
multivariate biorthogonal functions. A complete list of all possible limits is given in loc. cit. Here we just 
like to remark that this includes multivariate versions of all orthogonal polynomials in the (?-Askey scheme, 
and the correspondence of a-vectors to families of polynomials in the g-Askey scheme is given by the table 
in [Ij Section 7]. In the next section we consider another interesting case, the Pastro polynomials. 

6. Pastro Polynomials 

In this section we study the special points in the degeneration scheme for which the limits are families of 
biorthogonal polynomials (i.e. outside the q-Askey scheme where we have orthogonality). These polynomials 
are multivariate analogues of the polynomials studied by Pastro |8j . Specializing the parameters in the correct 
way gives us the only example of orthogonal polynomials on the unit circle (as opposed to the real line for 
the g-Askey scheme) contained in our degeneration scheme: the Macdonald polynomials [5]. 

The top level of these polynomials is llllpp, associated to the vector a = (— |, 0, |, i; 0, i) For the 

first function F6^\zp~i;top~'i ,ti,t2P^ ,t3pi;uo,uip^) the valuation turns out to be zero, and we get three 
different representations, corresponding to the expansion (^0)) with respect to the parameters v — to, v — t2, 
and V = ts, which make the evaluate by the bulk difference equation. Note that the expansion with v = ti 
does not give a nice limit, as in this case, the limiting interpolation function are inside the octahedron, and 
hence do not depend on z. Thus we get the following expansions (you can find the definitions of the limiting 
binomial coefficients and interpolation functions in Appendix [K[ in particular are just the Macdonald 
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polynomials) . 



Ef^\ '^M (^)^°(ti'^o) p*(«) / -1 , N «(M').-2n(M) ^ ,n-l, 



Clit^-'toh) \hUoJ E2,l^„_l^^^^ ] C°(t")C0(t"-1^3»l) 

Univariately these three expansions correspond to 



(if one considers the last expansion), and expressions of this series as 302, and as a 3(^)1, which are related to 
each other by H (III.6) and (111.7)]. 

Notice that the parameters of the function only appear in certain combinations. In particular if we write 



.1/2 



A = f_n-it-^uo ' ^ ~ t^-Hsui ^^'-^ ~ t'i-H2t3Uiz- '^^^ define (using the homogeneity of the Macdonald 
polynomials R^^l^^) 



' CO(i^)^ ^ ,lkW£;2jf] CO(t")CO(|) ^■'^^ B 

and note 

1/2 

p1")(z.; to : ti, ^2, ^3; «o, -1) - P^"^ (,„-4,3^^,^ ; ^;r4l^- F^^^' 

For the right hand family (the functions with uq and ui interchanged) we notice that by the symmetries 
from Lemma |3.2[ and the permutation symmetry in the t^'s we have 

rt(n) / _1 _1 11 1 \ " (n) / 111 _1 1\ 

R\'[zp i;top *,ti,t2P'^,t3p2;uip2,uo) = R\'(zp'i;top^,tip2,t2P S ^3; ""i, ""oP" ) 

n(ri)/ -1-1,1,1,-1, 1 

= R\'(z p 4;^op^^lp^^2P st3;ui,uop2) 

Rx'^ {z^^p-- i ; t2P~ Kh, topi , tipi ; m , ugp^ ) 

" A°,(P-^^ I ph-Hoh,ph--Hoh,p-'^j^,pi^) 

As a corollary we see that the limit on the right hand side is up to a constant equal to the limit on the left 
hand side with different parameters. Indeed we have 

11 111 / 1 \ 

Qx'^ ■■= lc(^R'}!^\zp^i;toP~Ktl,t2pKt3P^;Uipi,Uo)j = L2(«-l)^^^^^^^^ j P\'\zr^yt2,t3,to,ti]Ui,Uo), 

and the corresponding valuation is 0. In particular, setting 

qi"\wf,A,B) :=p[-\±;B,A), 

we have 



Q\ '{zf,to : ti,t2,t3;uo,ui} = q\ H ,„-i, , ^ T^TT;^ ' T^T^TT ■ 

Note that this is the same parameter correspondence as between P\ and p\. As already shown in T_ , the 
univariate instances and qf^ are equal to the Pastro polynomials (3.1)]. 

We would like to consider a special case. Indeed we want to specialize t^uiV^^^ 1, or equivalently 

(n) 

B ^ q. The simplest expression to do this in is in the expansion of in terms of the Macdonald 
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polynomials R^'^ , and to use px. We cannot just substitute B = q as for example the term C^{q/B) in the 
denominator would vanish. However, taking the limit p — )■ in |10[ (4.3)] with a — >■ ap" (and < a < 1) 
shows that 



= <5A,M(-a<z)l^|-|^l 

6=1 



so plugging this in we see that the sum reduces to a single term and we get 

(n) 

thus up to a constant these are just the Macdonald polynomials. If we also want Q\ to also be the 
Macdonald polynomials, we should also specialize uotit"''^ 1, or equivalently A ^ q. 
The limit of the biorthogonality relation gives us 

{p(-\.;U;ur;t,q),Q(;\-;U;ur;t,q))^S,,^t-'"^'\huo)-^'\^^^ 

where 



n\{t; q)^ j-J-^ {qP "■/tiUo,qP "/ia^^i; 7c" i<j<7:<„ 



X 



-A- 6l(gt^ "zj/totiUQ;q) dzj 



or in terms of px and g;!^ this becomes 

(pi")(.;Ai?),.l,")(-;Ai.)) 

where 

/r \ {Q;qT TJ^*^'^^^ [ ff \ t \ TT {wj/wk,Wk/wj;q) 



(t; g)« A A (AtJ-i, BtJ-i; (?) ^ i<.Ki<„ (^w^j'M, ^w^fcM; 9) 



l<j<fe<ri 

-Q 0{q^^'^Wj;q) dwj 



X 
3 



{Awj /q^/'^, B/wjq^/'^;q) 2'!riwj 



The contours here are chosen with conditions similar to the condition in Theorem 13. 3[ that is, C = C ^ is a 
deformation of the unit circle, which includes tC and the poles at Zj = tQq^^°, respectively Wj — Bq^2+^>«^ 
and excludes the poles at Zj — t2^q^-°, respectively Wj — A^^q^~^^i° . In the univariate case this inner 
product relation reduces indeed to the biorthogonality [8, (3.2)] of the Pastro polynomial^ Observe that 
the specialization A, B ^ q turns this measure into the measure of the Macdonald polynomials, in particular 
the univariate part of the measure is in that case just the constant 1. 

We can also take the limit in the evaluation duality ([22l) . There are two distinct evaluation dualities we 
can get in the limit, we can evaluate at a sequence based around to, or evaluate at a sequence based around t2 
(evaluating at ti or would not lead to an equation for Pastro polynomials) . Taking the leading coefficient 
in 

^^"^(<o<""''?'''p"^;top-3 : ti,t2pKt3P^;uo,uipi) = ^<f)(toi"~'q^'p~3;ioP"^ : i'l, ^2^^ , ; "o, UiP^) 
leads to 

)(tor~*g«' ;to,ti,t2,h;uo,ui) = Pj,"\iot"-'q^' ; to, ^1,^2,^3; uo, "1), 



'^Pastro has a necessary condition on the parameters which we lack, the difference is that he insists that the contour is the 
unit circle, while we look at more general contours (which are allowed to make detours to include points outside the unit circle 
and exclude points inside the unit circle) 
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where the dual parameters are given by 

X I T~ " totr ^ toUr 

to '■= ytotlt2t3/q, tr := Ur '■= — — . 

to to 

This is equivalent to 

At'"'~'^q'^* At''^~^q^* 

The equation based at t2 leads to 



This second evaluation duality, when specialized at B = q, reduces to the evaluation duality of the Macdonald 
polynomials. 

In this case it turns out that the evaluation duality gives an equation between the same functions (evaluated 
at different points). In general we would get an evaluation duality where both sides of the question are 
different functions (i.e. the dual parameters are not associated with the same limiting family of biorthogonal 
functions). 

Let us finally consider the difference operators. First we consider the equations = i?^"'' where 

the parameters of the two .R^"^'s are related by a g-shift. Expanding Dq"^ as a sum of 2" terms one quickly 
sees that the valuation of the individual terms only depend on cr e {±1}" only through \a\ = 
particular there are three options: Either the term with cr = (1,1,..., 1) dominates, or the term with 
a = (—1, —1, . . . , —1) dominates or all terms have the same valuation. In the latter case we have to worry 
about cancellation, where the valuation of the individual terms might be lower than the valuation of the 

(n) ~ (n) 

sum. It turns out this docs not happen, which we can easily check as we know the valuation of Dq R\ . 

For example, taking the difference equation for D'^\uo,to,ti) and taking the direct limit in each of the 2" 
terms leads to 

'■'■t—^ u v / o-e{±l}" i:o-i = l j:(7i=— 1 

= P^'^\zi;tQ : ti,t2,t3;UQ,Ui), 

/(|<j-|+n)/2\ 

where the t^ ^ > comes from the cross terms with Ui = <Tj = 1. Recall that there are 6 = (2) different 

difference operators like this, one for each pair of t-paramctcrs. The difference equations look simpler when 
expressed in the rcscaled polynomials px so in the list below we will use those. 



• D^''\uo,ta,ti): 



^^((M..)/.) ^ ^ l^:^p(")(^,,(3-.)/V/,,B,)=pl")(^,AB), 



1 ■ 1 - - WjIWi 
-lj:<7j=l 



• L>^"-'(uo,to,i2): 



(Te{±l}" J:cr,=l 



xiC'^'t--) n n ^^^Pi"^K.(^--'/v,55)=pi")(«.,A5) 
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• D^\uq , to , t^) , in this case the term with = 1 for all i is dominant in the limit, and the resulting 
"difference" equation of one term is the trivial equation p^"^ {wi ;A,B) = p^^^ {wi ; A, B) . 

• i'g"'' {uo,ti , ^2): In this case we have to be careful about the different normalization if we interchange 
to and one of the other t- variables. This accounts for the g"'^' term on the right hand side. 



;:(Ti = -l j:(Tj=l 



• D^\uQ,t2,t^). Again we have one dominating term, and the equation becomes trivial. 
Curiously enough, if we consider the first and fourth of these equations, we have two difference operators 
with different shifts of the parameter (either by respectively (g, 1)) acting on the same function 

p"^\-\ A/ Bq) giving the same result up to a generalized eigenvalue. This is in sharp contrast to the 
situation for general biorthogonal functions in our scheme, in which generalized eigenvalue problems only 
arise from products of pairs of such operators. Another remark is that the fifth difference equation preserves 
i3, so we can specialize B = q on both sides and obtain a difference equation for Macdonald polynomials. 
This operator is Macdonald's original difference operator (3.2)], of which the Macdonald functions are 
eigenfunctions. 

We have two more difference operators, namely D- and Z?+. These also have proper limits and we get as 
limit for that 

1 0"/2 



V _!_l_('_n(-«-kl)/2^-5(;)-(n-l)(k|+n)/2+(<l-l+"'/=^) 



11 /^. PA+i^g'^ '^'^/'wf,A/q,B)=\[ ^ p[>iw.;A,B) 



The operator moreover gives us 

n V TT (1- ^ ) TT (1-^) 



,,(<M.")/^) ^ i^^p(")(,(--)/V;A,s,).pi';),„(-.;AS). 

l:(7i — — \j:(Tj—\ 

Appendix A. Limits of interpolation functions 

In what follows we are going to define c^^x^k for different values of e. These are the coefficients used in the 
branching rules for different limits of the interpolation functions R* ^ . This convention allows us to define 
these limiting interpolation functions simultaneously. 

Definition A.l. For any e for which c^,x,k is defined we define the corresponding interpolation function 
R*''^\zi; a,b; q,t) recursively by setting -R*^o^ — 1 and for A 7^ setting -R**\^ = 0, and using the branching 
rule 

(23) Rl^+'\...,z,,...,v;a,b;q,t)^ ^ c,,x,.R:^:\. ■ ■ , z,, . . . ;a,b;q,t). 



As an important part in the definition of Ce^A,K tlie following expression often arises. 
Definition A. 2. For a y^' X we set (^), , = 0. For n ^' \ we aenne 

^t^' [a.,t]-q,t ^ ■' 

The definition of this new binomial coefficient was inspired by the following lemma. 
Lemma A. 3. Let /.i X. Then we find for < a < 1 

val((^) UO, val((^) ] =a{\X\-\^,\), 

^^,/A\ ^_ Cl{t)C+{a) (-i^)l^lq"(^')i-2n(A) y;^ 



Proof. It is convenient to first rewrite the elliptic binomial coefficient as 
A\ _ C^{t)C+{a) -j-j eiq'^^-H^'.-'pq^q^^^-H^+^^'.-'-p) 



^'l [aA,<i,t,v C^it)Cna/t) ^^^f-^^ e{q^^^-H^^'-^pq,q>^^-H'+^'^-^-p) 

X'.—it' 

U{iJK,.diq^'^'-'t^'^''''pqci,q^^+'-'t^-'''^''a;p) 

KM 

Y\[^,J)ex^{q^'^+'-^t-^''^-'pqa,q^^+^-H^-^'^~'a■p) 
KM 

~ {~a)M + t2«(A) + |A| C^(t)C+{a/t) ^} 0(qM.-J + ltM;-i^qA,-jil+A;.-^.p) 

A' — u' 

X ^^^^^ 

n(.,j)GA0(g^-+^"i-''^-*a,g^'+^-ii'-^^-^a;p) 

The first identity is obtained by plugging in the definitions ([2]), respectively ([S]), of and as products 
over boxes in A, and distributing those products over the terms 11(1 j)eA A' =/i' ^'^'^ Y\{i ^p,' ■ should 

be observed that G A | A^- = /x^} = {(i, G | A^- = ^^ j. The second equality then follows by 

using 9{px;p) = —^d{x;p) for all theta functions which visibly contain a p in the argument. The resulting 
prefactors can be simplified to the term given by some combinatorial arguments. The two trickiest ones are 
that 

(i,j)GAi i l<j<Pi i (i,j)eA 

and 

E ^^-^ E (^^-1)= E E E (^-1)- E {^-^) = n{\)-n{^,). 

In the latter calculation we use that if A^ 7^ /i^ then A^ ^ fi'j + 1 (as /i ^' A), and that if £ X/fi then 
Aj = j as must be the lowest box in its column (again as /i -<' A). 
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We have now written the elhptic as a product where the leading coefhcients of the terms are clear, so 
the two limits follow immediately. □ 

The case (^)jgj is closely related to the coefficient in the branching rule for Macdonald polynomials, indeed 

Lemma A. 4. We have 

\Ai/[o] 

where ipx/fi as given in [SI page 341] oj Macdonald' s book. 

Proof. Notice that in the notation of [5] we have 6^(s) — b\{s) whenever s G A — Rx/^ — Ca//j- Thus we can 
rewrite ipx/^i in terms of our notation as 



^a(s) ^.-"-i, bx{s) 1 



V ^a(s) &a(s) \Ai/roi' 



□ 



This normalization is very useful in order to prove functions defined by the proper branching rules are 
monic. 

Proposition A. 5. Suppose clonic, x.k = (k)[o]^"(^' where is a monic polynomial in v of degree 

A„, then the coefficient of — z^^ • • • z^" in the resulting interpolation functions R* (i^) monic x ^I'^O'^^ ^• 

Proof. The proof hinges on two observations. The first is that there is exactly one chain of partitions, which 
leads to terms z^, so the coefficient of z^ can be explicitly written as 



i),(Al 



The second is that if A = /i • fc, i.e., Xi = /i^ for i < n and A„ = k (hence k < /i„_i), then (^)jj^j = 1- Indeed 
in this case we see that if {i,j) G A with A^- = /i^, then i < n, so A^ = /i^. Thus it follows that all terms in 
the defining product of (^)jpj equal 1. □ 

In particular this proposition allows us have a well-defined normalization for all polynomial interpolation 
functions given below: we normalize them so that they are monic. 

Now we are ready to explicitly write down all limits. As mentioned in Section |4] the limits under consid- 
eration are associated to the faces of different dimension of the tessellation with tiles as in Figure [TJ Below 
we list all those faces (up to symmetry) and give the limiting interpolation functions. 

A.l. The vertex. All the vertices give essentially the same limit by the shifting formulas. For simplicity 
we use the vertex (0, 0, 0). 

Definition A. 6. We define 

C°(raz;±i)C°(^^^i) ^l^l 



^"'''''^ \^/[*"fl;«,*C'0(i"az;±i)(70(f«±i) 

The definition of cv^x.k automatically provides us with the interpolation functions Ry^ and we get the 
following proposition. 

Proposition A. 7. We obtain the limit 

val{Rl^"'\zi; a,b; q,t; p)) = 0, 

lciRl^-\zf,a,b;q,t;p))=R<l\z,,aMq,t f^^^^^^ 
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Proof. The valuation val{c\^i^{a,h;v)) — and the value of the leadmg coefficient lc{cx^f^{a,b;v)) follow 
directly from the expressions for the binomial coefficient fLemma lA.3[) and of the A^'s (which follow from 
(llip ). The limit for the interpolation function then follows from the branching rule using induction to n. 

It could be that the newly defined Ry^ is identically zero, in which case the valuation of the interpolation 
functions would be more than 0, the value given here, and the expression for the leading coefficient would 
be completely wrong. Fortunately we have Corollary lA. 151 below, which shows this is not the case. □ 

The proofs for the other limits of the interpolation functions are essentially the same as this one, so we 
will omit them from now on. To get an idea of the kind of coefficients we have here, it should be observed 

that 

n (1 

'^^ ' (■ij)eA/K 

so it is a polynomial in v of degree |A| — with leading coefficient 

So in particular C'°(i"au*^)/C'°(i"au*^) is a polynomial in + of degree |A| — For the terms involving 
b we notice 

thus we get a rational function with a polynomial in w + of degree Ai + |k| — |A| in the numerator and 
the given polynomial (f w^^; q)xi of degree Ai in the denominator. 

A. 2. Edges. For the edges we use the names El, E2, E3 and E4, which correspond to 

El The edge connecting (0, 0, 0) and 5, ^); 
E2 The edge connecting (0,0,0) and (^,—5, ^); 
E3 The edge connecting (0,0,0) and (0, 1,0); 
E4 The edge connecting (0, 0, 0) and (1, 0, 0). 

The edges E3 and E4 are symmetric under the change ( — ^ — C, while El and E2 change into the final two 
edges (mod lattice translations). 

Definition A. 8. We define the coefficients 

^/A\ g^(t"g)C°(^) ,|.| 



CE3,X,k — 



A\ (^^(rawii) (j-"(^')r(^)(-ar)-l^l 
K/ [0] C'0(t"aw±i) g-"(«')i"W(-a<»)-l''l 



Proposition A. 9. We have for < a < 1 for the valuations 

val{Rl^''\z,p''/^;ap'^^^,bp''/^;q,t;p)^0, val{R*^''\z,p'''^-ap"'^ ,bp-'''^-q,t-p) = a\\l 

val{Rl^''\z,p'/^;ap^^^,bp'/^-'';q,t;p) ^ a\X\, m/(i?*^"^(z,; ap", fe; g, t;p) = 0, 



and for the leading coefficients 

Note R*^2\ does not depend on a, so we left it out of the notation. Similarly R*^^\ is independent of 

b. For R*^i \ and R*^2 \ slightly more subtle. These satisfy an invariance of the following form: for 
arbitrary x e C* we have 

(24) R*^'^^^{ziX] ax, bx; q, t) = a, 6; q, t), R*^^^^{ziX; ax, b/x; q, t) = R*^^^^{zi; a, b; q, t). 

Thus, in essence, they too have one less variable than Ry^- 

We would also like to remark that we normalized R*^^\ to be a monic polynomial in Zi + z^^ . The 

interpolation function R*^^\ is moreover equal to Okounkov's |B] i3C„-symmetric interpolation function, via 

R*E^\i'^i'^ ~ {o,t^^^)^^^ P\{zi / at^ ; q, t, a), as can be easily seen by comparing the branching rule given here 
with 16, (5.2)]. 

A. 3. Triangles. The different triangles are called Fl, F2, F3, and FA (for face). 

Fl The triangle with vertices (i,— i,i), (0,0,0) and (5,5,^); 
F2 The triangle with vertices (0,0,0), (0,1,0) and (i,i,-i); 
F3 The triangle with vertices (5,-5,5), (0,0,0) and (1,0,0) ; 
F4 The triangle with vertices (0,0,0), (1,0,0) and (5,5,-5); 

Note we get all eight different triangles (up to lattice shifts) when we reflect these four by C — > — C- 



Definition A. 10. We define the coefficients 



CF2,\,K = 



CF3,\,k 



C_F4,A,K — 



A\ Cl{t"av) (-ra)-l^lg-"(^')t"(^) 
K/ CO{t^av) (-t"o)-l«lg-"('^')t"('=) 

•A\ '^°(bl7) ^,|K|-lA| 



Proposition A. 11. We have for < a and < /3 and a + [3 < 1 

val{Rl^"\z,p"/^+^/^; ap^/^+P/^, 6p"/2-/3/2. ^. ^ ^j^j yal{Rl^"\z,p-^/^; ap''/^ bp^/^+l"; q, t; p)) = 
^;a^(i?I'"^(V/2;ap"+^/^6p-^/2;g,^;p)) =/3|A|, m?(i?*(") (z.p-^/^; ap"+^/^ 6//^; g, i;p)) = 
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2\ l^'l 



lc{R<-\z.p-^'^; ap-^^'\ 6//^; f^p)) = Rt\{z^; b; q, t)^^ (-^) 9"^"'^^-^"^"' 



Like for Rei and Re2, these functions all satisfy a shifting equation, so they have one less variable than 
apparent. These equations are 

Rp^ ^ {xZi , (2x) ^ Rp-^ y^i^Zi^ (l)^ R p2,X {"^Zi ^ (xj — ^ ^F2 A ' ^) ' 

R*F3^xi^Zi;b/x) = R*j!:!^\{zi;b), R*^2\{^zf,b/x) = x^^^'^ R*j2\izi;b). 

The interpolation function R*pi \ correspond to the interpolation functions studied by Okounkov in [7]. 
The relation is R*^'^\{z,; a) = P{{z,/fa){V'~^a)~\^\. 

A. 4. Tetrahedra. We denote the tetrahedron with vertices (0, 0, 0), (1, 0, 0), (5,^,5) and (i, — i, i) by T. 
All tetrahedrons are shifts of this tetrahedron, or shifts of the reflection of this tetrahedron in the plane 

C = o. 

Definition A. 12. We define the coefficients 

\'*/ [0] 

Proposition A. 13. We have for < a, /3, 7 and a + /3 + 7 < 1 

These functions associated to the interiors of the tetrahedra are the famous Macdonald polynomials [5]. 
Proposition A. 14. Let P\ denote the Macdonald polynomials. Then we have 

Px{zi;q,t) ^ R'^^\zi;q,t). 

Proof. The Macdonald polynomials can also be defined using a branching rule, indeed using 'W, (7.9)' and 
(7.14)'] we see that the coefficients in the branching rule for the Macdonald polynomials are exactly ipx/fj,- 
By Lemma lA.41 those are the coefficients in the branching rule for Rt,x- The result now follows using 
induction. □ 

As a corollary we obtain the following important result. 

Corollary A. 15. The families of interpolation functions {^**'^''}AeP„ for e equal to one of V, El, E2, E3, 
E4-, Fl, F2, F3, F4, T form a basis for the symmetric functions in Zi. 

Proof. The family {R*r^^x}xePn forms such a family as they are identical to Macdonald's interpolation func- 
tions; the other families all degenerate to this family by taking the right limits, so they must also form such 
bases. □ 
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A. 5. Octahedron. The interior of the octahedron consists of two square pyramids and the square separating 
these pyramids. This gives us three hmits, however all of them are independent of z^. As functions in Zi they 
are thus not interesting, however we can give a explicit expression for the limiting constant, which means we 
obtain some nontrivial combinatorial equations. 
Let us label the different polytopes by 

S The square with vertices (0, 0, 0), (^, ^, ^), (1, 1, 0) and (^, ^, —\); 
PI The pyramid with base S and apex (0, 1, 0); 
P2 The pyramid with base S and apex (1, 0, 0). 

Definition A. 16. We define the coefficients 

/A ,-\-\ /^\ /^\ i-M 



N"-/ [t"f] \'^/ [0] \'*/ [0] 

Proposition A. 17. We have the following evaluations 

i?9xfa b-a t)^ ^n(A)^-(n-l)|A| 

^A ^A y^) 

Moreover we have for a,/?, 7 with \a\ + \/3\ + I7I < 1 the limit 

valiRl^''^ {z,p^/'; ^^1/2+0/2+7/2^ ^^1/2+0/2-7/2. ^. ^ 

Proof. First we obtain an expression for the limits of the interpolation functions using the branching rule as 
before. The case 7 = corresponds to the square, while 7 < corresponds to PI and 7 > to P2. This 
gives us the valuations of the interpolation functions in these cases, and expressions for the leading coefficient 
in, respectively, Rs.x, Rpi,x and Rp2.\- In particular we see that the leading coefficient is z-independent. 
But we also know the evaluation (Proposition 12. 6| ) for general interpolation functions 

Rt\xt--^-aM - A",(^ I t--'ax, ^). 

Plugging in a — >■ ap^/^"'""/^"'"^/^, b 5pi/2+a/2-7/2^ g^j^^j 2; xp^l"^ in here, shows that the valuation on the 
right hand side is and the leading coefficient 1. Thus we get another expression for the leading coefficient 
for special values of Zi. However, as the leading coefficient is (by inspection) independent of the z^, it must 
be equal to this constant. The resulting equality gives us the evaluations of Rs.x, Rpi,\ and Rp2.\. D 

The value Rpi^x is equal to a special value of the Macdonald polynomials. Indeed, given a chain = 
A(o) ^' ^' ... ^' = A, the coefficient of Hi ^1^'''^^''"" ' in the Macdonald polynomial Pa equals 
rii (a('»-i))[o]' Rpi this product of binomial coefficients appears with a factor 



J|t|A*'-'I^J|tEr..|A*'=VA''=-'|^-Q, 



Thus we obtain that Rpi^x — Pt.aC^"^^, ^"^^i • • ■ : 1), and its evaluation thus corresponds to the famous 
principal evaluation of Macdonald polynomials. Likewise the evaluation of Rp2,x also corresponds to this 
same evaluation of Macdonald polynomials. 

A. 6. Binomial coefficients. Recall that the binomial coefficients were defined in terms of the interpo- 
lation functions. Obtaining the limits of the binomial coefficients is simply a question of plugging in the 
corresponding limits of the interpolation functions. These limits were all already discussed in |10[ Section 



Proposition 12.81 (which says the binomial coefficients are p-elliptic in a and b) tells us that we only need 
consider limits for for a,/3 modulo the lattice 1? . Notice that the binomial coefficients only 

[ap" .^p'^] ;g,i;p 

use the interpolation functions with a = C,, that is, in the plane at the front of Figure[TJ The relevant picture 
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(0,1) = (0,1,0) V 




v(i,i) = (i,i,i) 



(0,0) = (0,0,0) V E2 v(l,0) = (i,-i,i) 



Figure 2. The different limits for the binomial coefficients. Points are labeled by {a, (3) — 
(a/2, /3- a/2, a/2) 

of the different limits is given in Figure [H The second labeling is included as it corresponds to the one in 
Figure [T] notice that we only have the front of the parallelepiped and reflected it left-right. 

Definition A. 18. We define 



A 

V/^/ El,[x];q,t 

\^^/ E2,[b]-q.t 
v/^/ S3,[a];q,t 
vA^y Fl-q.t 
vM/ F2:q,t 

for n > /(A), /(^). 



r'O/' l^ 



C^(q)Cl{aq) 



R*^^\{q>^H^-\t^-) 



^F2.uy- 



Proposition A. 19. We have the following limits for < a < /3 < 1 



lim 



P^O \tJ.J [a,b];qA;p \f-J V,[a,b];q,t 



lim 



P^O [ap'',bp"];qA;p \l^J El,[a/b];q,t 



lim 



p->0 [ap'',b];q,t;p E2,[b];q,t 



lim 



lim 



P^O [ap'3,bp°]:g,t:p V^/ Fl:q,t 



P^O [a,6p°];g,t;p E3,la]-q,t 

A\ /A^ 



lim 

P^'^ \l^J [ap" ,bpl^];q,t:p \l^JF2:q,t 



Proof. We plug in the limits of this section in Definition 12.71 We simplify the results a bit (by using 
invariances of the form (|24|) ) to indicate that, for example, {^) j ^j^j ^ only depends on a/b (as opposed to 
a and 6). □ 
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